MPI-PhT/94-43 

TUM-T31-65/94 

June 1994 



Parity Violating Longitudinal Muon Polarization 

in K + — > 7r + /i + /i~ 

o^ ■ Beyond Leading Logarithms Q 

ov 

^ • Gerhard Buchalla and Andrzej J. Buras 

o: 

CN ■ Max-Planck-Institut fur Physik 

t— i ' - Werner-Heisenberg-Institut - 

CNj , Fohringer Ring 6, D-80805 Miinchen, Germany 

m ■ 

o. 

^r | D-85748 Garching, Germany 

o ,' Abstract 

We generalize the existing analyses of the parity violating muon polarization 
asymmetry A^ in K + — > 7i + fi + u.~ beyond the leading logarithmic approxima- 
tion. The inclusion of next-to-leading QCD corrections reduces the residual 
dependence on the renormalization scales, which is quite pronounced in the lead- 
ing order. This leads to a considerably improved accuracy in the perturbative 
calculation of the short distance dominated quantity Alr- Accordingly this will 
also allow to obtain better constraints on the Wolfenstein parameter q from fu- 
ture measurements of A LR . For -0.25 < q < 0.25, V cb = 0.040 ± 0.004 and 
m t = (170 ± 20)GeV we find 3.0 • 10~ 3 < \A LR \ < 9.6 • 10~ 3 . 



Technische Universitat Miinchen, Physik Department 



* Supported by the German Bundesministerium fur Forschung und Technologie under con- 
tract 06 TM 732 and by the CEC science project SC1-CT91-0729. 



It has been pointed out by Savage and Wise JTJ that measurements of muon po- 
larization in K + — > n + fi + fi~ decay can give valuable information on the weak mix- 
ing angles and in particular on the parameter g in the Wolfenstein parametrization 
[0]. Indeed as shown in |j], [|, the parity- violating asymmetry 

I ^lr | = | ^^ |= r- | Re£ | (1) 

ifi + U 

is dominated by the short distance contributions of Z-penguin and W-box dia- 
grams with internal charm and top quark exchanges, while the total rate is com- 
pletely determined by the one-photon exchange amplitude. The interference of 
this leading amplitude with the small short distance piece is the source of the 
asymmetry Alr- Here Tr (Tl) is the rate to produce right- (left-) handed /x + , 
that is /i + with spin along (opposite to) its three-momentum direction. The fac- 
tor r arises from phase space integrations. It depends only on the particle masses 
rriK, Tn-K and m M , on the form factors of the matrix element (tt + | (sd)y-A I K + ), 
as well as on the form factor of the K + — > 7r + 7* transition, relevant for the 
one-photon amplitude. In addition r depends on a possible cut which may be 
imposed on 9, the angle between the three-momenta of the fi~ and the pion in 
the rest frame of the fi + fi~ pair. Without any cuts one has r = 2.3 ||. If cos# 
is restricted to lie in the region —0.5 < cos 9 < 1.0, this factor is increased to 
r = 4.1. As discussed in ||, such a cut in cos 9 could be useful since it enhances 
Alr by 80% with only a 22% decrease in the total number of events. 
Re£ is a purely short distance function depending only on CKM parameters, the 
QCD scale A.jjg and the quark masses m t and m c . We will discuss it in detail 
below. 

A-lr as given in (|1|) has also been considered by Belanger et al. ||], who 
emphasized its close relation to the short distance part of the decay amplitude 
Kl — > fi + fi~. Unfortunately the authors of ref. @ did not include the internal 
charm contribution to A^. As we will show explicitly below the charm contri- 
bution cannot be neglected as its presence increases the extracted value of g by 
roughly A^» = 0.2. 

Let us briefly summarize the theoretical situation of Alr- The "kinemati- 
cal" factor r can be essentially obtained from experimental input on the particle 
masses and the form factors. The form factor / describing the K + — > 7r + 7* vertex 
has been discussed in detail in H within chiral perturbation theory. While the 



imaginary part can be reliably predicted, the real part is only determined up to a 
constant to be fitted from experiment. On the other hand, data on K + — > ir + e + e~ 
decay || allow to extract the absolute value of this form factor directly. We will 
adopt this approach, following ||. Since the imaginary part of / is quite small 
[0, M, we then also have the real part Re/. In principle Im/ could yield an extra 
contribution in ([[]) proportional to Im£. We have checked, based on the approach 
of || that this contribution is below 1% of the dominant part shown in ([]]) and 
can therefore be safely neglected. Clearly the factor r involves some uncertainty 
due to the experimental errors in the form factors, which can however be further 
reduced by future improved measurements. For the present discussion we will 
assume fixed numerical values for r. 

Besides the short distance part of A LR there are also potential long distance con- 
tributions coming from two-photon exchanges, which have also been discussed by 
the authors of j3]. These are difficult to calculate in a reliable manner, but the 
estimates given in || indicate that this contribution is substantially smaller than 
the short distance part, although it cannot be fully neglected. Therefore the short 
distance effects are expected to safely dominate the quantity Alr and we shall 
concentrate our discussion on this part, keeping in mind the possible existence of 
non-negligible long distance corrections. 

The short distance physics leading to A LR is generally considered to be very 
clean, as it can be treated within a perturbative framework. However this does 
not mean that it is free of theoretical uncertainties. An indication of the involved 
error due to the necessary truncation of the perturbation series in the strong 
coupling constant a s can be obtained by studying the sensitivity of a physical 
quantity to the relevant renormalization scales on which it should not depend 
in principle. The existing short distance calculations of Alr |j], ]3|, [| include 
QCD corrections in the leading logarithmic approximation (LLA) p|. As it turns 
out, they suffer from sizable theoretical uncertainties due to the residual scale 
dependences. The main purpose of our letter is to extend the analyses of |], [|. [§ 
beyond the leading logarithmic approximation thereby reducing considerably the 
theoretical uncertainties in question. To this end we will use our next-to-leading 
order analysis of K L — ► /i + /i~ presented in [g, 0]. 

Our discussion of the Cabibbo-Kobayashi-Maskawa matrix will be based on 
the standard parametrization [Hj, which can equivalently be rewritten in terms 



of the Wolfenstein parameters (A, A, g, 77) through the definitions [TT1 

s 12 = A s 23 = AX 2 Sl3e -^ = A\ 3 (Q-iri) (2) 

The unitarity structure of the CKM matrix is conventionally represented through 
the unitarity triangle in the complex plane with coordinates (0, 0), (1, 0) and (g, 
fj) where 



g + irj = 



V udV* b 



v cd v; b 

To better than 0.1% accuracy g = g(l — A 2 /2) and f] = 77(1 — A 2 /2). 

Following || [| it is straightforward to generalize the expression for Re£ of 
beyond leading logarithms. We find 



(3) 



Re£ = k 
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(4) 



(5) 



Here A = \V US \ = 0.22, sin 2 9^ = 0.23, x t = m 2 /M^ and X { = V*V id . The 
function Y, relevant for the top contribution, is given by 



O'.s 



Y(x) = Y Q (x) + -rY^x) 

4:71 
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where x M = /i 2 /M^ with fi = \i t = 0(m t ) and 
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The QCD correction Y\ has been calculated in M. Next 

Y NL 



Pn 



A J 



(9) 



(10) 



where Ynl represents the renormalization group expression for the charm contri- 
bution in next-to- leading logarithmic approximation (NLLA) calculated in 0. It 
reads 

Ynl — Cnl — B NL (11) 

where Cnl is the Z-penguin part and B NL is the box contribution, relevant for 
the case of final state leptons with weak isospin T 3 = —1/2. We have 



x[m) „|i 

NL 



C NL = -4^A7 5 
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where 



a,(M w ) a s (//) , , 

A = r~r- K c = — 7— r (loj 

a s (/i) a s (m) 
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A+ = A^s K_ = K^r R 33 = K^ (14) 



(-1/2) _ x(m) || 

n NL — 7 ^c 



*-*> terw^ 



, u 2 329 15212 _ 30581 _ 

In -Hr 1 A 2 H A A 2 

m? 12 625 7500 



(15) 



Here A 2 = A _1//25 , m = m c , x = m 2 /Myy- In ( |i"2|) - ( [TBI) the two-loop expression 
has to be used for a s (/z) and /i = /x c = 0(m c ). The explicit /i-dependences in 
the next-to- leading order terms (||), (|T2|) and ([T^) cancel the scale ambiguity of 
the leading contributions to the considered order in a s . The consequences of this 
feature will be discussed later on. Numerical values of Pq are given in table 1 
where m c = m c (m c ). 

It is evident from (|]) and (f|) that, given |A ii? |, one can extract ReA t : 



\*lr\ (l-^ 

Re\ t = -A 5 ^ A 2 ^ " (16) 



i-£m 



Y(x t ) 

Furthermore, using the standard parametrization of the CKM matrix we obtain 
from the definition @ 



1 + 4si 2 Ci 2 ReA t /s| 3 - (2s 12 c 12 Im\ t /s?, 3 ) 2 - 1 + 2s 2 . 
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2c 2 s 2 (1?) 

ZC93S1 



-23*12 





P 


A I7s> m c [GeV] 


1.25 


1.30 


1.35 


0.20 


0.132 


0.141 


0.150 


0.25 


0.135 


0.145 


0.154 


0.30 


0.139 


0.148 


0.158 


0.35 


0.142 


0.152 


0.162 



Table 1: The function Pq for various A^j^ and m c . 

Up to the very accurate approximations that V C( iV* b is real (error below 0.1%) 
and C13 = 1 (error less than 10~ 5 ) (|T7|) is an exact relation. Using the excellent 
approximation ImA t = 7]A 2 \ 5 [O, we see that a measured value of ReA t deter- 
mines by means of ( |TTD a curve in the (g, ?7)-plane. Since the dependence on 
ImAj is however very small, this curve will be almost parallel to the 7/-axis. Thus, 
knowledge of | Alr\, hence ReA i; implies a value for g (or g) almost independently 
of ImA t . For simplicity we shall neglect ImAt in ([17]) completely, which introduces 
a change in g of at most 0.01. It is evident that the more general treatment can 
be easily restored if desired. Note that the charm sector contributes to g the 
non-negligible portion 



A^ charm « P /{A 2 Y{x t )) « 0.2 



'M 



In order to demonstrate briefly the phenomenological consequences of the 
next-to-leading order calculation we consider the following scenario. We assume 
that the asymmetry A^# is known to within ±10% 



A LiJ = (6.0 ± 0.6) ■ 10 



-3 



(19) 



where a cut on cos#, —0.5 < cos# < 1.0, is understood. Next we take (rrii = 
rhi{mi)) 

m t = (170±5)GeU m c = (1.30 ± 0.05)GeU V cb = 0.040 ± 0.001 (20) 

A^ = (0.30 ± 0.05)GeU (21) 

The errors quoted here seem quite reasonable if one keeps in mind that it will 
take at least ten years to achieve the accuracy assumed in (|i~9|). The value of mt 







A(A LiJ ) 


A(m t ) 


A(V cb ) 


A(m c ) 


A(A^) 


Q 


-0.06 


±0.13 


±0.05 


±0.06 


±0.01 


±0.00 



Table 2: £> determined from A^# for the scenario described in the text together 
with the uncertainties related to various input parameters. 

in ( p0|) is in the ball park of the most recent results of the CDF collaboration 

& 

In table 2 we have displayed the central value for g as it is extracted from A^ 

(flTEp and (|I7|)) in our example, along with the uncertainties due to the parameters 

involved. This is intended to indicate the sensitivity of g to the relevant input. 

The combined errors due to a simultaneous variation of several parameters may 

be obtained to a good approximation by simply adding the errors from table 

2. It is interesting to compare these numbers with the renormalization scale 

ambiguities, which inevitably limit the accuracy of the short- distance calculation. 

If we vary the scale in the charm- and in the top sector as lGeV < fi c < 3GeV 

and 100GeV^ < fi t < 300GeV, respectively, keeping all other parameters at their 

central values, we obtain the following range for g 



-0.15 < g< -0.03 (NLLA) 

-0.31 <g< 0.02 (LLA) 
We would like to emphasize the following points: 



(22) 
(23) 



The error in g from (0), which illustrates the theoretical uncertainty of 
the short distance piece alone, is not negligible. It seems however moderate 
when compared to the errors shown in table 2. We stress that (^) is based 
on the complete next-to-leading order result for Re£. If only the leading 
log approximation is used instead, the range obtained for g is by almost a 
factor of 3 larger (|23"D . 



The error in (p2|) is almost entirely due to the charm sector. Indeed, if we 
vary only yU c , keeping fi t = m t fixed, the corresponding interval for g reads 
(—0.14, —0.03). This illustrates once more, that the charm sector, being 
the dominant source of theoretical error in the short distance contribution 
to Alr, should not be neglected. 

6 



• In the case x<l, which is relevant for the charm contribution, the function 
Y has a very special structure. Expanding the renormalization group result 
Ynl to first order in a s one finds (here x = rn^/M^r) 

Y NL = ^ + £-x In 2 x (24) 

I 47T 

We observe that the leading logarithms ~ x In x, present in the Z-penguin- 
and the box part, have canceled in Y NL , leaving the subleading term x/2 as 
the only contribution in the limit a s = 0. On the other hand QCD effects 
generate an a s x\n 2 x "correction", which is of the order O(xlnx), hence 
a leading logarithmic term! As a first consequence the charm function Y 
is enhanced considerably (by a factor of ~ 2.5) through strong interaction 
corrections, compared to the non-QCD result. (This feature is in a sense 
similar to the case of the rare decay b — > sj.) A second point is that the 
x/2 term, though formally subleading, is important numerically. Working 
within LLA one is then faced with the problem of how to deal with this 
term since it should strictly speaking be omitted in this approximation. Let 
us illustrate this issue in terms of the g determination in our above example. 
We find g = —0.12 if we use the LLA formulae (with fi c = m c , fi t = rn t ) 
and simply add the x/2 piece. By contrast, omitting this term and using 
the strict leading log result we obtain g = —0.20. The scale ambiguities 
are very similar in both cases, roughly three times as big as in the next-to- 
leading order discussed above. For definiteness we have included the x/2 
part to obtain (|23|) . 

The problem of the x/2 term is naturally removed in the next-to-leading 
logarithmic approximation (Yjvl) where this contribution is consistently 
taken into account. 

Finally we give the standard model expectation for A^, based on the short 
distance contribution in (|j), for the Wolfenstein parameter g in the range —0.25 < 
g < 0.25, Vc = 0.040±0.004 and m t = (170±20)GeV. Including the uncertainties 
due to m c , Ajjg, /i c and /i t and imposing the cut —0.5 < cos9 < 1, we find 

3.0 • 10~ 3 < \A LR \ < 9.6 • 10~ 3 (25) 

employing next-to-leading order formulae. Anticipating improvements in V c b, 
m t and g we also consider a future scenario in which g = 0.00 ± 0.02, Vcb — 

7 



0.040 ± 0.001 and m t = (170 ± 5)GeV. The very precise determination of q used 
here should be achieved through measuring CP asymmetries in B decays in the 



LHC era [p^| . Then fl25|) reduces to 

4.8 • 10" 3 < |A M | < 6.6 • 10~ 3 (26) 

In both of the scenarios the lower (upper) limit for A^ would be smaller by 
0.6 • 10~ 3 (1.3 • 10~ 3 ) if the charm contribution was omitted. 

In this letter we have generalized the short distance calculation of the muon 
polarization asymmetry Alr in the decay K + — > 7r + /i + /! - to next-to-leading 
order in QCD. We furthermore discussed the theoretical uncertainties involved 
in this analysis. We have demonstrated that the complete next-to- leading order 
calculation achieves a reduction of the rather large scale ambiguities in leading 
order by a factor of ~ 3 and is necessary to provide a satisfactory treatment of 
A LR . This is particularly important since long distance contributions to A LR 
seem to be small, though perhaps not fully negligible [§]. 

In any case a measurement of Al R would yield a very interesting and useful piece 
of information on short distance flavordynamics and the unitarity triangle which 
is worth pursuing in future experiments. 
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